First detection of glycolaldehyde outside the Galactic Center 

M. T. Beltran 2 , C. Codella 3 , S. Viti 4 , R. Neri 5 , R. Cesaroni 6 

ABSTRACT 

Glycolaldehyde is the simplest of the monosaccharide sugars and is directly 
linked to the origin of life. We report on the detection of glycolaldehyde 
(CH 2 OHCHO) towards the hot molecular core G31.41+0.31 through IRAM PdBI 
observations at 1.4, 2.1, and 2.9 mm. The CH 2 OHCHO emission comes from the 
hottest (> 300 K) and densest (>2xl0 8 cm" 3 ) region closest (< 10 4 AU) to 
the (proto) stars. The comparison of data with gas-grain chemical models of 
hot cores suggests for G31.41+0.31 an age of a few 10 5 yr. We also show that 
only small amounts of CO need to be processed on grains in order for existing 
hot core gas-grain chemical models to reproduce the observed column densities of 
glycolaldehyde, making surface reactions the most feasible route to its formation. 

Subject headings: ISM: individual objects: G31.41+0.31 — ISM: molecules - 
stars: formation 



1. Introduction 

Glycolaldehyde (CH 2 OHCHO), an isomer of both methyl formate (HCOOCH3) and 
acetic acid (CH3COOH), is the simplest of the monosaccharide sugars. The importance of 
this organic molecule is that it can react with propenal to form ribose, a central constituent of 
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RNA, and is, therefore, directly linked to the origin of life. Hence, it is crucial to investigate 
the occurrence of glycolaldehyde in the universe, especially in star-forming regions, where 
stars associated with planetary systems are expected to form. Interst ellar glycola l dehyd e 
has first been detected towar ds the Galactic cent er so urce Sgr B2(N) bvlHollis et al.l (120001 ) 



and successively confirmed by lHalfen et al.l (120061 ) and iReauena- Torres et al. 



J2008h . through 



single-dish observations. Interferometric observations (iHollis et al.ll200ll ) of the Sgr B2(N) 
region, known as the Large Molecule Heimat (LMH), have shown that, unlike its isomers 
methyl formate and acetic acid, glycolaldehyde has a large spatial scale, 60". 



Beltran et al.l (120051 ) have reported the detection of a bright line at 220.46 GHz, observed 



with the IRAM Plateau de Bure Interferometer (PdBI), towards the massive star- forming 
regions G31.41+0.31 and G24.78+0.08. In the latter case, such a line has been detected 
in two distinct cores, Al and A2. This emission has been interpreted as possibly arising 
from the transition 202,i8~193,i7 of CH2OHCHO. However, it should be noted that a single 
tentative dete ction of an interste llar molecule transition is not enough for establishing its 
identification (jSnyder et al.ll2005l ). We have performed an observational campaign to detect 
multiple transitions of CH2OHCHO in G31.41+0.31 in order to confirm for the first time the 
presence of glycolaldehyde in a compact core in a star forming region, possibly associated 
with the material assembling the newly born massive protostars. 

G31.41+0.31 is a very massiye star -forming; region located at a distance of 7.9 kpc 
( jChurchwell. Walmsley. fc Cesaronil Il990l ). associated with a hot (T ~ 300 K) molecular 
core (HMC), with a lumino s ity of ~3xl0 5 L P) and thought to be heated by one or more O-B 
(proto) stars. IBeltran et al.l (120041 ) have detect ed a clear v e locity gradient perpendicular to 
the direction of a bipolar outflow imaged by lOlmi et al.l (119961 ) at the center of the core 
(see Fig. [2]). Such a velocity gradient has been interpreted as rotation of a massive toroid 
(M gas ~490 M & ) about the axis of the outflow. In addition, recent SMA interferometric 
observations at 226 GHz have detected CN and H 2 CO in red-shifted absorption towards 
the dust continuum of the hot core (Beltran et al., in preparation), which indicates the 
presence of accretion in such a massive toroid. Therefore, G31.41+0.31 very much resembles 
the 20 M Q young stellar object G24.78+0.08 Al, where infall, outflow, and rotation have 
been simultaneously detected for the first time towards a massive (proto)star (IBeltran et al. 
2006h . 



In this letter, we report the first detection of glycolaldehyde towards a HMC located 
outside the Galactic center, and show that there is a viable mechanism for glycolaldehyde 
to be abundant in HMCs by means of chemical modelling. 
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2. Observations 

We carried out observations at 1.4, 2.1, and 2.9 mm with the PdBI on March 2003, 
June 2008, and July 2008, respectively. Line data have a spectral resolution of 3.4, 0.33, and 
0.45 km s _1 at 1.4, 2.1, and 2.9 mm, respectively. Channel maps were created with natural 
weighting, obtaining synthesized beams and la RMS of l'.'lxO'.'S and 50 mJy/beam/channel 
at 1.4 mm, 4'. / l x3'.'2 and 13 mJy/beam/channel at 2.1 mm, and 5 / .'4x4'. / 4 and 9 mJy/beam/channel 
at 2.9 mm. 



3. Detection of glycolaldehyde 



Figure 1 reports the spectra observed at 1.4, 2.1, and 2.9 mm towards the peak position 
of the HMC G31.41+0.31. As can be seen in the top panel, besides the CH 3 CN (12 8 -11 8 ) 
and CH-3 13 CN (126-lle) transitions, other lines are detected. In particular, there is a clear 
emission peak at 220465.86 MHz. Although the line is blended with the CH 3 CN (12 8 -11 8 ) 
and CH 3 13 CN (12 6 -11 6 ) transitions, there is no doubt of the presence of a separate line 
at 220465.86 MHz. On the basis of the Jet Propulsion Laboratory, Cologne, and Lovas 
databases for molecular spectroscopy, the only possible identification is the CH 2 OHCHO 
(20 2 ,i 8 -19 3 ,i7) line at 220463.87 MHz (£ u =120 K; S/j, 2 = 65.39 D 2 ). The difference between 
the observed and the laboratory frequencies (~1.99 MHz), albeit non negligible, is less than 
the spectral resolution of 2.50 MHz. Such a discrepancy is not surprisin g since it has been 



found in several g lycolaldehyde transitions observed towards Sgr B2(N) (iHollis et al.l 12000 



Halfen et al.ll2006l ) and could reflect the uncertainties on the CH 2 OHCHO rest frequencies 
in the spectral catalogues. As mentioned above, the glycolal dehyde 20^8-19-^7 line is 
blended with two molecular transitions of CH3CN and CH3 13 CN. iBeltran et al.l (120051 ) fitted 
the three lines with Gaussians having the same line width (FWHM), and separations in 
frequency identical to the laboratory values (allowing for the above mentioned difference of 
1.99 MHz for the glycolaldehyde transition). Here, we prefer instead to leave all parameters 
free to allow for possible differences in the FWHM of different molecules and for the sake of 
consistency with the method adopted to fit the other two glycolaldehyde lines (see below). 
The new value of the FWHM of the 1.4 mm line is 9. 2 ±0.7 km s -1 . 

To confirm the tentative detection of glycolaldehyde, we checked for other CH 2 OHCHO 
lines that have already been detected towards the Galacti c center. Ba s ed on the Galactic 
center CH 2 OHCHO survey performed at 2 mm and 3 mm by lHalfen et al.l (120061 ). we observed 
towards G31.41+0.31 the CH 2 OHCHO (10i ; g-9 2i8 ) and (14 0il4 -13i ; i3) transitions respectively 
at 103667.907 MHz (E u =32 K; Sfi 2 = 26.22 D 2 ) and 143640.936 MHz (£ u =53 K; S/j? = 67.56 
D 2 ). These two transitions are the perfect follow-up for this search, because their excitation 
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is lower than that of the 220 GHz line and their rest frequencies are well known. In addition, 
the CH 2 OHCHO (14 0i i4-13i i3) is less affected by line blending with respect to the 1.4 mm 
one, whereas the (lO^g-92,8) hne is not contaminated by other spectral emissions. As can be 
seen in Fig. [TJ both transitions have been clearly detected at the expected frequencies. The 
derived FWHMs are 4.6±0.6 km s" 1 and 7.1±0.2 at 103.67 and 143.64 GHz, respectively. 
Both values are significantly less than the width of the 1.4 mm lin e (9.2±0.7 km s" 1 ) . 



FWHMs spanning of up to a factor of 2 have also been observed by lHalfen et al.l (120061 ) 
towards the Galactic center. In our case the fact that the line width increases with the 
energy of the transition could be due to an excitation effect, if the higher energy lines arise 
from a smaller, more turbulent region. We can hence conclude that we have obtained the 
first detection of glycolaldehyde outside the Galactic center. 

The integrated intensity of the 1.4 mm line, J T-q dv = 201.8±5.9 K km s _1 , is definitely 
higher with respect to those of the 2.1 mm (4.8±0.1 K km s _1 ) and 2.9 mm (1.0±0.1 K 
km s _1 ) lines. This is the consequence of the beam dilution effect. Figure 2 shows the maps 
of the three transitions of CH2OHCHO inte grated emission as well as the velocity gradient 
detected in CH3CN by IBeltran et al.l (120041 ) towards the center of the hot core. As seen in 
the maps, the emission of glycolaldehyde comes from the central region, where the millimeter 
continuum peaks. Our angular resolution at 1.4 mm (l'.'lxO'.'S) is sufficient to barely resolve 
the CH2OHCHO emission (bottom- left panel of Fig. 2), thus minimising beam dilution 
effects. On the opposite, the angular resolution of the 2.1 and 2.9 mm maps, ~ 4"-5", causes 
a significant beam dilution, thus reducing the measured brightness temperature. 



Unlike what has been observed towards the Galactic center (IHollis et al.ll200ll ). the 
CH2OHCHO emission does not have a large spatial scale but is rather concentrated towards 
the center of the core. In fact, the glycolaldehyde emission has a deconvolved size at the 
50% of the peak of the emission of ~1'.'3 (~10300 AU), which is de finitely smaller tha n 
that derived from standard HMC tracers, such as e.g. methyl cyanide (IBeltran et al.ll2005l ). 
which is of ~2'.'4 (~19000 AU). As a c onsequence, as show n by the rotational temperature 
and column density maps (see Fig. 7 of IBeltran et al.ll2005l ). CH2OHCHO traces the hottest 
(T > 300 K) and densest (tih 2 >2x10 8 cm -3 ) gas in the surroundings of the newly formed 
star(s), which is located closer to the peak of the continuum emission. 

Given the difference of the angular resolution of the three CH2OHCHO maps, the com- 
parison of the different lines has been done by convolving the 1.4 mm and 2.1 mm observation s 
to the 2.9 mm HPBW (5'.'4x4 / .'4). By using the expressions reported by IHollis et al.l (120001 ). 
we derived the rotational temperature diagram, finding that the diagram is quite flat, which 
does not allow to derive the column density from the fit. This is probably due to the high 
opacity of the CH2OHCHO lines. In fact, on the one hand, the brightness temperature 
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of the 1.4 mm line is quite high (~30 K). On the other hand, after smoothing all maps 
to the resolution of the 2.9 mm line, the brightness temperature of the three transitions 
turn out to be very similar. This suggests that all three lines are optically thick. There- 
fore the beam-averaged column densities obtained from the 2.9 and 2.1 mm lines assuming 
T=300 K (iVcH 2 OHCHO — 3xl0 15 cm -2 ), as well as that derived from the better resolved 
1.4 mm transition (~10 17 cm -2 ) should be considered as lower limits. A rough estimate of 
the CH 2 OHCHO abundance has been obtained by using the colu mn density determined 
from the continuum dust emission at 1.4 mm (IBeltran et al.l 120041 ) . The estimated abun- 
dance is of the order of 10~ 8 . Such a high value is consistent with the abundances found by 
Requena- Torres et al. r|2008h towards the Galactic center. 



Besides the detection of glycolaldehyde, we wish to demonstrate not only that its column 
density is consistent with the predictions of theoretical models, but also that other species 
involved in the same formation route are detected with column densities compatible with 
those predicted by the models. One of these is HCOOCH 3 . It appears that the HCOOCH3 
line at 220445.79 MHz cannot be used to estimate the column density because th e value 



obtain ed would be implausibly large (A^hcooch 3 ~10 cm ). In fact, according to iTurner 



( 119911 ). transitions of complex molecules associated with small (< 1 D 2 ) Sfi 2 , often exhibit 
intensities in excess by orders of magnitude with respect to the LTE values based on the 
large Sp? lines. We t hus decided to estimate the HCOOCH3 column density using data from 
Cesaroni et al.l (119941 ) for a transition with S/j 2 = 13 D 2 . The result is listed in Table HI after 
correcting the HCOOCH3 estimates for the differen t beam with respect to the C H^OHCHO 
observations. For this purpose, the expressions of iRequena- Torres et al.l (120061 ) have been 
used. Following the same approach, the column density of methanol (CH 3 OH; Table HJ, 
another molecule in yolved in the formation process of glycolaldehyde, has also been collected 
from the literature (iGibb. Wvrowski. fc Mundvll2003h . Note that if HCOOCH 3 and CH 3 OH 
were optically thick, the column densities could be much higher. In the following, we compare 
the column densities thus estimated (see Table [Q with those obtained from theoretical 
models. 



4. The origin of glycolaldehyde in HMCs 

In this section we briefly investigate whether the column densities of glycolaldehyde 
that we derive are indeed chemically feasible. Our purpose is to investigate whether current 
chemical models of HMCs can account for the observed abundance of glycolaldehyde in 
massive star forming regions. 

To date, several gas-phase and solid-phase routes of formation for glycolaldehyde have 
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reaction of HCO with methanol or with formaldehyde (or, more likely, via intermediate, e.g. 
CH 2 OH + HCO, where the first species is a direct product of methanol). While such reactions 
would be too slow in the gas phase to produce detectable abundances of glycolaldehyde in 
the HMCs lifetime, the high densities (10 7 cm" 3 ) of such region s may lead to fast surface 



reactions: it has been proposed for example ICharnley fc Rodgerd (120051 ) that such reactions 



may occur in close proximity via the hot secondary electron generated by the passage of a 
cosmic ray through the ice, or via photoprocessing of grain mantles by UV starlight which 
create a high concentration of radicals in the bulk interior of mantles. Grain-grain collisions 
then provide excess heat causing radical-radical rea c tions to occur and form large organic 
molecules (Sorrell 2001). In fact, iBennet fc Kaiser! (120071 ) have performed an experiment 



where methanol and CO ices were irradiated with energetic electrons to mimic the presence 
of cosmic ray ionizations and found that glycolaldehyde was formed provided that methanol 
and carbon monoxide are in close proximity in the ices. 

For the purpose of this investigation, we will adopt as route of formation of glyco- 
laldehyde surface reactions of HCO, H2CO and CH3OH, but we note that these somewhat 
arbitrary choices do not imply that other routes are not important. Also note that in this 
simple model, the formation of glycolaldehyde on ices via reactions of HCO with methanol 



i mplie s intermediate passages (just like the formation of methanol from CO in IViti et al. 



( 120041 ) which we assume to be fast. As far as the destruction of glycolaldehyde is concerned, 
gas-phase photodissociation would be highly inefficient in such environments due to the high 
visual extinctions. Destruction by surface reactions will be taken into account by varying 
the efficiency of formation. Glycolaldehyde formation via neutral-neutral reactions on the 
grains will affect the abundances of other HMC species such as CH3OH, and HCOOCH3. 



We have used a published HMC model (IViti et al.ll2004l ) modified to inclu de the forma- 
tion re actions for glycolaldehyde as well as other surface reactions summarized in lBottinelli et al 
( 120071 ) to form complex surface molecules (Table [2]). The model already included hydro- 
genation on grains. We have run a small grid of models where we varied the formation rate 
coefficients of the reactions forming CH 2 OHCHO and related species. Our best fit model 
is one where the efficiencies of the new reactions listed in Table |2] need not be large (e.g. 
only ~ 2% of CO needs to be converted in H 2 CO). The range of ages indicated by the best 
fit models is consistent with a typical age for a HMC (~10 5 -5xl0 5 yr). In absence of any 
experimental measurements, we do not know whether the efficiencies employed for Reactions 
2 and 7 (Table [2]) are realistic. Nevertheless, we note that the best fitting models imply rea- 
sonably small CO conversion efficiencies and probabilities of formation, probably feasible on 
surfaces where grains act as catalyst. More specifically both reactions occur among products 
of hydrogenation of CO hence the necessary mobility is low as all the radicals are already in 
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close vicinity to each other. 

5. Conclusions 

We report the first detection of glycolaldehyde towards a star forming region, the HMC 
G31.41+0.31 that is hosting massive young stellar object(s). G31.41+0.31, together with 
G24.78+0.08 Al, is the only massive core where evidence of infall, rotation, and outflow 
have been simultaneously detected. The maps of CH 2 OHCHO show that the emission, 
which is smaller than that derived from standard HMC tracers, is clearly associated with 
the most central part of the core. This emission traces the hottest (> 300 K) and densest 
(>2xl0 8 cm" 3 ) gas in the surroundings of the embedded high-mass protostar(s). Therefore 
we propose glycolaldehyde as an excellent tracer of the innermost regions of HMCs. To 
detect glycolaldehyde towards other HMCs it is crucial to carry out high-angular resolution 
observations because the interferometer will filter out the extended emission of other lines 
that might overlap with glycolaldehyde and affect its detectability. 

The detection of glycolaldehyde has also profound implications for the chemistry of star 
forming regions as it can help to constrain the evolutionary stage of the core, and therefore of 
the embedded O-B (proto) stars. In fact, we show that, for existing HMC gas-grain chemical 
models, the age of the HCM has to be a few 10 5 yr to reproduce the observed column densities 
of glycolaldehyde (as well as methyl formate and methanol). In addition, only small amounts 
of CO (~10-15%) need to be processed on grains, making surface reactions the most feasible 
route to the formation of glycolaldehyde. 

We thank Prof D. A. Williams for helpful discussion as well as a thorough reading of 
the manuscript. SV acknowledges financial support from an individual PPARC Advanced 
Fellowship. 
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Table 1: Source-averaged column densities towards the G31.41+0.31 HMC 
CH 2 OHCHO a HCOOCH 3 a CH 3 OH b 
(cm -2 ) (cm -2 ) (cm -2 ) 

>lxl0 17c 3.4xl0 18d ' e 7.5xl0 17c 



a Assuming a kinetic temperature of 300 K, as measured bv lBeltran et al.l (|2005f ). 
fc From lGibb. Wvrowski. fc Mundvl |2003h . 
c From the CH 2 OHCHO (20 2 ,i 8 -19 3 ,i7) data. 

d From the HCOOCH 3 -E (9 6 , 3 -8 6 ,2) (110652.89 MHz; E u =51 K) data of ICesaroni et alJ (|l994h . 
e In case of optically thick emission, the column densities could be higher. 
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Table 2: List of reactions added to the Viti et al. (2004) model. "M" denotes species in the 
solid phase. 



N 


Reaction 


1 


CO + 4(MH) =>> MCH 3 OH 


2 


CO + MCH3OH =>- MHCOOCH3 


3 


H 2 CO + MH MCH3O 


4 


MCH3O + MHCO MHCOOCH3 


5 


CO + 2(MH)^ MH 2 CO 


6 


CO + MH ^ MHCO 


7 


MH 2 CO + MHCO + MH =>- MCH 2 OHCHO 
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Fig. I. — Beam-averaged spectra in T B scale of the CH 2 OHCHO (20 2 ,i 8 -19 3 ,i 7 ), (14 ,i4- 
131J3), and (10i i9 -9 2 , 8 ) at 220463.87, 143640.94, and 103667.91 MHz, respectively, as ob- 
served towards the central position of the G31.41+0.31 hot core. Rest frequencies are 
pointed out by vertical bars. Upper panel (from Beltran et al. 2005): The glycolalde- 
hyde line is blended with the CH 3 CN (12-11; K=8) line. Two additional lines are present: 
(i) 13 CH 3 CN (12 6 -11 6 ; labeled by K'), and (ii) HCOOCH 3 -A (2511,15-269,18) (2204445.79 
MHz; i? u =272 K) which could contain an emission contribution due to the CH 2 OHCHO 
(18 4 ,i4-17 4 ,i3) (220433.51 MHz; £ u =108 K) line. The continuous line shows the fit to the 
group of three lines formed by the CH 2 OHCHO (20 2 ,i 8 -19 3 ,i 7 ), CH 3 CN (12-11; K=8) and 
13 CH 3 CN(12-11; K'=6); dotted lines draw the three individual Gaussian curves used for the 
fit. Middle panel: The CH 2 OHCHO line is part of a spectral pattern containing also the 
HCCC 13 CCN (54-53) (143636.63 MHz; £ u =183 K), C 2 H 3 CN (33 2 , 3 i-32 4 , 28 (143646.50 MHz; 
£ u =620 K), and C 2 H 5 OH (29 2 , 28 -28 3 , 26 (143651.78 MHz; £ u =415 K) lines. The results of 
the fit as drawn as in the Upper panel. Lower panel: Besides the glycolaldehyde emission, 
an unidentified spectral pattern is present around 103674 MHz. Solid curve shows the fit of 
the isolated CH 2 OHCHO line. 
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Fig. 2.— Map of the intensity integrated under the CH 2 OHCHO (10 ,9-9 2 ,8) line at 103.67 
GHz (a), the CH 2 0HCH0 (14 0il4 -13 lil3 ) line at 143.64 GHz (6), and the CH 2 0HCH0 
(20 2 ,i8-19 3)1 7) line at 220.46 GHz (c), towards the hot molecular core G31.41+0.31. The 
3 a contour of the CH3CN (12-11) emission averaged under the K=0, 1, and 2 components 
is shown in light blue, (d) Close-up o f the central region that shows the velocity field of 
the toroid in G31.41+0.31 mapped by iBeltran et al.l (120041 ) . The cross marks the position 
of the 1.4 mm radio continuum source, and the dots the 7 mm radio continuum sources 



(lAraya et al.ll2008l ). The synthesized beam is shown in the bottom left. The contour levels 
are: (a) from 0.01 to 0.09 in steps of 0.02 Jybeam -1 ; (b) from 0.04 to 0.52 in steps of 0.08 
Jybeam -1 ; and (c) from 0.10 to 0.94 in steps of 0.12 Jybeam -1 . 



